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The family of binary compounds including TaAs, TaP, NbAs, and NbP was recently discovered 
as the first realization of Weyl semimetals. In order to develop a comprehensive description of the 
charge carriers in these Weyl semimetals, we performed a detailed and systematic electronic band 
structure calculations which reveal the nature of Fermi surfaces and their complex interconnectivity 
in TaAs, TaP, NbAs, and NbP. Our work report the first comparative and comprehensive study 
of Fermi surface topology and band structure details of all known members of the Weyl semimetal 
family and provide the necessary building blocks for advancing our understanding of their unique 
topologically protected low-energy Weyl fermion physics. 


I. INTRODUCTION 

The discovery of Dirac fermions as low-energy quasi¬ 
particle excitations in graphene and on the surfaces 
of topological insulators has drawn significant atten¬ 
tion in both fundamental physics research and device 
applications^ — . Since these Dirac-like fermions prop¬ 
agate as massless relativistic particles they behave dif¬ 
ferently from the conventional charge carriers in met¬ 
als, semiconductors, and insulators. Recently, a new 
form of massless fermions with lifted degeneracy at the 
nodal point, from four-fold to two-fold, has been pro¬ 
posed to exist in condensed matter systems through 
a time reserval or inversion symmetry breaking mech¬ 
anism. This symmetry breaking operation modifies 
the quasi-particle’s dispersion relation from a Dirac to 
a Weyl equation^ - — Many theoretical proposals ex¬ 
ist for realizing Weyl semimetals that possess inter¬ 
esting physical properties, such as discontinuous Fermi 
arcs and negative magnetoresistance due to the chiral 
anomalyA ^ " 10 ! 19 " 29 ! 31 ! 32 ! 34 A Weyl node with definite chi¬ 
rality is associated with the Berry curvature and may be 
thought of as realizations of magnetic monopoles in mo¬ 
mentum space j 3 Q i 61 

Recent theoretical works have proposed the realiza¬ 
tion of the Weyl semimetal state in the inversion sym¬ 
metry breaking TaAs family ] n i 13 Shortly after the pre¬ 
diction, the first Weyl semimetal was experimentally dis¬ 
covered in TaAaiS. The electronic Weyl semimetal state 
in TaAs was experimentally observed through photoe¬ 
mission spectroscop y 12 ! 14 . Electrical transport experi¬ 
ments have shown that TaAs has very high mobility^ 


and, therefore, is consistent with the protected nature of 
the Weyl fermions and reported signatures of the chiral 
anomal y 36 ! 37 . Soon after the initial experimental discov¬ 
ery other independent photoemission experimental works 
have confirmed the Weyl semimetal state in TaAs and, 
furthermore, the Weyl state in other members of the same 
family, which includes NbAs and TaP^^r— 

Provided the inclusion of spin-orbit coupling, these 
four non-centrosymmetric compounds are proposed to re¬ 
alize Weyl fermions without breaking time-reversal sym¬ 
metry and, consequently, are easily studied under ambi¬ 
ent conditions. For example, no external magnetic field 
or pressure is required for the realization of the Weyl 
semimetal state in this TaAs family ] 11 ! 12 Conventional 
semimetals are materials that have a small overlap in 
energy between the valence and conduction bands. By 
contrast, a Weyl semimetal has valence and conduction 
bands that touch at a set of discrete points in the bulk 
Brillouin zone, called the Weyl nodes. In the presence of 
additional doping, Fermi pockets, or multiple Weyl nodes 
at different energies, a Weyl semimetal will naturally con¬ 
tain electron and hole charge carriers, giving rise to a 
compensated semimetal. This is indeed the case for the 
TaAs family of Weyl materials as reported in photoemis¬ 
sion and transport experiment s 12 i 14 ! 35 ~ — The Fermi sur¬ 
faces can originate from the conduction or valence band 
crossing the Fermi energy, the Weyl cones, and their com¬ 
binations. The charge carriers are expected to have small 
concentrations with a high mobility, at least for the Weyl 
fermions. The interplay between tunable concentrations 
of multiple charge carriers, high Fermi velocities, and the 
protected nature of the Weyl fermion carriers suggests a 
new avenue for designing novel electronic devices. The 
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Weyl semimetals are also good candidates for exploring 
critical phenomena near the Fermi energy, due to their 
small carrier concentrations and linear band touchings. 
Therefore, it is of interest to explore the electronic struc¬ 
ture of these four compounds to better understand their 
similarities and differences. 

Previous theoretical studies of this family focused on 
establishing TaAs as a Weyl semimetal and therefore only 
reported a few key aspects of the electronic structure for 
that purpose. On the other hand, many systematic de¬ 
tails of the band structure, which are crucial for spec¬ 
troscopy and transport studies, have not been studied. 
Furthermore, a comparative study of the electronic struc¬ 
ture of TaAs with the other members remains lacking. In 
this study, we perform a comprehensive first-principles 
calculation to explore the properties of Fermi surfaces for 
the entire family of Weyl semimetals, TaAs, TaP, NbAs, 
and NbP, and, thereby, providing the necessary building 
block for future investigations of Weyl semimetals from 
a theoretical, experimental and application perspective. 
The paper is organized in the following manner: In Sec¬ 
tion HD we provide the computational details. In Sec¬ 
tion mu electronic band structures and density of states 
are studied to improve our basic understanding of these 
four compounds. In Section UYl we discuss the details 
of the Fermi velocities of Weyl fermions and Fermi sur¬ 
faces. Finally, in Section [V] we present a brief discussion 
covering the important conclusions reached. 


II. COMPUTATIONAL DETAIL 

First-principles calculations of TaAs, TaP, NbAs, and 
NbP, which have space group Mimd (109), a body- 
centered-tetragonal structure, were performed using the 
OpenMX code. The OpenMX code was based on 
norm-conserving pseudopotentials generated with multi 
reference energies^ and optimized pseudoatomic basis 
functionS i 53 i 54 For each Ta atom, three, two, two, and 
one optimized radial functions were allocated for the s-, 
p-, d-, and /-orbitals (s3p2d2/l), respectively. For As, 
Nb, and P atoms, s3p3d3f2 , s3p3d3fl , and s3p3d2fl 
were adopted, respectively. A cut-off radius of 7 Bohr 
was adopted for Ta, Nb, and P basis functions while a 
9 Bohr cut-off radius was used for the As basis func¬ 
tions. The spin-orbit coupling was taken into account 
through j-dependent pseudopotentials ^ and the gener¬ 
alized gradient approximation (GGA) was adopted for 
the exchange-correlation energy functional ] 56 i 57 The cut¬ 
off energy of 1000 Ry was used for numerical integra¬ 
tions and for the solution of the Poisson equation. For 
the quality of /c-point sampling, a 17 x 17 x 5 mesh for 
the conventional unit cell was adopted. The experimen¬ 
tal lattice parameters were chosen in the calculations and 
are listed in Table HI The conventional unit cell and first 
Brillouin zone of the primitive unit cell of the represen¬ 
tative TaAs is plotted in Fig. [lj 

Throughout this paper, for brevity, we will label the 


TABLE I. Experimental lattice parameters (c/. Fig. [T]). 


TaAs 

TaP 

NbAs 

NbP 

a (A) 3.437 

3.318 

3.452 

3.334 

c (A) 11.656 11.363 

11.679 

11.376 

u/c 0.333 

0.334 

0.333 

0.334 



FIG. 1. (Color online), (a) Conventional unit cell of TaAs. 
(b) Primitive first Brillouin zone of TaAs. TaP, NbAs, and 
NbP share the same body-centered-tetragonal structure with 
different lattice parameters listed in Tabled] 

crystal momentum in units of (27r/a, 27r/a, 2tt/c) unless 
otherwise specified. 

III. ELECTRONIC STRUCTURE 

In this section, the electronic band structure and den¬ 
sity of states for the four Weyl semimetal compounds 
are presented. With the combination of Ta, or Nb, and 
As, or P, very similar electronic structures can be found 
and, therefore, illustrates a flexibility in realizing Weyl 
fermions among these four elements. 

A. Symmetry 

The relevant symmetry elements in the TaAs family 
(space group IAimd) is the fourfold screw rotations along 
z -axis and two mirror reflections with respect to x- and 
y- axis. Due to lack of inversion symmetry, the rotation 
axis cannot sit perpendicular to the mirror planes. The 
combination of screw rotations and mirror reflections re¬ 
sults in diagonal glide planes. Because the screw rotation 
and glide reflection are compound symmetry operations 
of translation and point group operations, for specific sur¬ 
face termination the surface states will not contain these 
compound symmetries. For the bulk, however, due to 
translation symmetry, the energy spectrum displays the 
C 4 and mirror symmetries. In k space, there are mir¬ 
ror planes at k x = 0, k y =0 and k x db k y = 0. The k x —0 
and k y =0 mirror planes are important for the realiza¬ 
tion of Weyl nodes. On the mirror plane, each band 
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can possess a mirror eigenvalue of either +1 or —1 (+i 
or —i) in the absence (presence) of spin-orbit coupling, 
respectively. Bands with opposite mirror eigenvalues are 
allowed to cross, and those with the same mirror eigenval¬ 
ues are only permitted to anti-cross. The other symmetry 
is time-reversal, which is not a spatial symmetry. Since 
time-reversal symmetry relates k and —k states, its com¬ 
bination with the twofold rotation results in the energy 
spectrum being symmetric with respect to the k z = 0,1 
planes. 

Since a Weyl node is characterized by a chirality of 
d=l, which is defined by the integration of Berry curva¬ 
ture, a pseudovector, over a closed k -space surface (in 
unit of 27r) j 3Q i 61 , Weyl nodes related by mirror symme¬ 
try will possess an opposite chirality while those related 
by rotation and time-reversal symmetries will have the 
same chirality. Thus, a mirror pair of Weyl nodes, as 
mirror images, will possess, relative to each other, an 
opposite chirality, and two Weyl nodes at equal k x and 
k y , but opposite k z , will have the same chirality. Fur¬ 
thermore, the closed surface can be extended to a Fermi 
surface, yielding a Chern number for the Fermi surface^ 
and, therefore, reflecting the net chirality of the enclosed 
Weyl nodes. 



B. Band structure 

The electronic band structures of TaAs, TaP, NbAs, 
and NbP with and without spin-orbit coupling are shown 
in Figs.[2j[3j|U and[5j respectively. As expected, the band 
structures of these four Weyl semimetal compounds are 
characterized by the d orbitals of Ta or Nb atom, and the 
p orbitals of As or P atom are similar to one another near 
the Fermi energy. In the absence of spin-orbit coupling, 
the valence and conduction bands cross and form closed 
rings that are bounded on the mirror planes (k x = 0 and 
k y = 0 planes), due to the fact that both bands having 
opposite mirror eigenvalues. For the four Weyl semimetal 
compounds the nodal rings on the k x = 0 mirror planes 
are shown in Fig. [ 6 ] Similarly, the location of nodal rings 
on the k y = 0 mirror planes can be found by applying 
a C 4 rotation. Again, note that away from the mirror 
planes, which lack the necessary symmetry requirements, 
accidental band crossings are not found. 

In the presence of spin-orbit coupling, the aforemen¬ 
tioned nodal rings no longer exist. The valence and con¬ 
duction bands become fully gapped along the high sym¬ 
metry lines, as shown in Figs. [2j [3j SJ and [5j The gap 
opening is smaller when replacing the Ta atom for an 
Nb atom, which indicates weaker spin-orbit coupling in 
NbAs and NbP and is consistent with the strong spin- 
orbit coupling in heavy elements. 

The disappearance of nodal rings follows with the ap¬ 
pearance of Weyl nodes. Pairs of Weyl nodes around the 
rings with opposite chirality on the two sides of the mir¬ 
ror planes are found i 11 ^ 13 For each ring in Fig. [ 6 j three 
pairs of Weyl nodes are generated, and their locations, 


FIG. 2. Electronic band structure of TaAs (a) without and 
(b) with spin-orbit coupling. 



FIG. 3. Electronic band structure of TaP (a) without and (b) 
with spin-orbit coupling. 
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FIG. 4. Electronic band structure of NbAs (a) without and 
(b) with spin-orbit coupling. 



FIG. 5. Electronic band structure of NbP (a) without and 
(b) with spin-orbit coupling. 



FIG. 6. (Color online). Line nodes on the k x —0 plane formed 
by the crossing of valence and conduction bands in (a) TaAs, 
(b) TaP, (c) NbAs, and (d) NbP without spin-orbit coupling. 
The red solid circles indicated the projection of Weyl nodes 
on the k x =0 plane after turning on the spin-orbit coupling. 
For each circle, two Weyl nodes of opposite chiralities on two 
sides of the k x —0 plane are found. 


projected on the k x = 0 plane, are indicated by the solid 
circles. The location of the other pairs of Weyl nodes can 
be obtained by applying a C 4 rotation, which uncovers 
24 Weyl nodes for the first Brillouin zone. The detailed 
k -space coordinates and energies of two representative 
Weyl nodes are given in Table HU We denote Weyl nodes 
by W\ for those at k z = n (multiples of ^), and by W 2 
otherwise. The energies of W\ are lower than those of 
W 2 for all of the studied compounds. For a Weyl node 
below (above) the Fermi energy, an electron (hole)- like 
Fermi surface will enclose a Weyl node, respectively. 

Although the nodal rings are destroyed by spin-orbit 
coupling, new nodal rings are allowed to be created since 
the closest valence and conduction bands possess oppo¬ 
site mirror eigenvalues, =b i. In other words, bands can 
cross on the plane with suitable hopping parameters. In 
NbAs, such a ring is found at k ~(0, 0.41, 1.47) with a 
tiny radius < 0.005. The negligible k -space area is in 
accordance with the energy overlap of only < 0.2 meV for 
the ring. The energy of the ring is located approximately 
5.4 meV above the Fermi energy. Since the energy over¬ 
lap is tiny and occurs above the Fermi energy, we neglect 
the significance of this overlap in this study. The fact 
that spin-orbit coupling can annihilate or keep the nodal 
lines indicates that these lines depend on the presence 
of mirror symmetry but are not protected by it. How¬ 
ever, parity of the number of nodal lines is conserved 
because of topology. Since the nodal line is in a Z 2 class 
(1 + 1 = 0 )^ interaction between two nodal lines can 
produce zero or two lines. 
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TABLE II. The coordinates (in units of reciprocal lattice vec¬ 
tors of conventional unit cell) and energies (in units of eV) 
of two representative distinct Weyl nodes denoted as Wi and 
W 2 are given below. In each compound the energy of W 2 is 
higher than that of Wi . =b stands for a mirror pair of Weyl 
nodes. 



coordinate of 

Wi 

energy of Wi 

TaAs 

(±0.0072, 0.4827, 

1.0000) 

-0.0221 

TaP 

(±0.0074, 0.4809, 

1.0000) 

-0.0531 

NbAs 

(±0.0025, 0.5116, 

1.0000) 

-0.0322 

NbP 

(±0.0028, 0.5099, 

1.0000) 

-0.0534 


coordinate of 

w 2 

energy of W2 

TaAs 

(±0.0185, 0.2831, 

0.6000) 

-0.0089 

TaP 

(±0.0156, 0.2743, 

0.5958) 

0.0196 

NbAs 

(±0.0062, 0.2800, 

0.5816) 

0.0042 

NbP 

(±0.0049, 0.2703, 

0.5750) 

0.0259 


C. Density of states 

TaAs, TaP, NbAs, and NbP have similar electron con¬ 
figurations in the outer shells, namely d 3 s 2 and p 3 for 
Ta, or Nb, and As, or P, respectively. In order to fill the 
outer-most p orbitals of As and P, As 3- , P 3_ , Ta 3+ , and 
Nb 3 + are expected according to the ionic-bonding pic¬ 
ture, which should leave the Ta d or Nb d orbitals having 
a higher energy than the full-filled As p or P p orbitals. 
In Fig. [7J the partial density of states for Ta d, Nb d, As 
p, and P p are presented. One can directly observe that 
the d orbitals are strongly hybridized with the p orbitals. 
Although p character leads to the major contributions 
in the bottom band, the p orbitals are not fully occu¬ 
pied. In fact, the occupation numbers obtained in the 
basis of linear combination of atomic orbitals give only 
approximately half-filled p orbitals and the p orbitals also 
contribute to the energy above the Fermi energy. There¬ 
fore, it is inadequate to neglect the effect of p orbitals 
near the Fermi energy. 

A simple picture to understand the obtained density of 
states is that of the Ta, or Nb, d orbitals hopping strongly 
to the As, or P, p orbitals and giving rise to the bonding 
and anti-bonding bands across the Fermi energy. The 
Ta, or Nb, s orbital also participates in the hybridiza¬ 
tion and does not donate itself completely. Having this 
hybridization scheme, the semimetal feature is formed, 
which shows the valence and conduction bands separated 
by a valley shaped density of states at the Fermi energy. 
The density of states are not zero at the Fermi energy 
due to the lack of particle-hole symmetry and the two 
distinct Weyl nodes (Wi, W 2 ) possessing different ener¬ 
gies for each compound. Furthermore, the valence and 
conduction bands are found to cut through the Fermi 
energy and, therefore, forming hole and electron pockets, 
respectively. Note that the major contribution to the 
density of states around the Fermi energy is of d char¬ 
acter. The bonding and anti-bonding p orbitals are split 
into much lower and higher energy. 

Some additional information can be obtained from 



FIG. 7. (Color online). The orbital contributions of Ta d 
or Nb d and As p or P p in the density of states with spin- 
orbit coupling for (a) TaAs, (b) TaP, (c) NbAs, and (d) NbP 
are presented. The d and p contributions are colored red 
and green, respectively. The total formula contributions are 
plotted by black lines. 


Fig. H TaAs and TaP share a wider bandwidth rela¬ 
tive to the narrower bandwidth of NbAs and NbP, which 
reflects that 5 d electrons are more delocalized than 4 d 
electrons and that the replacement of As p by P p has 
less of an effect on the overall density of states. However, 
it is clearly observed that the density of states of P p in 
TaP is slightly lower than As p in TaAs. In addition, 
the density of states of P p in NbP is lower than As p 
in NbAs. This indicates that P p orbitals lose more elec¬ 
trons and are further away from the picture of fully-filled 
p orbitals. Another observation is that Nb d orbitals in 
NbP show the largest area below the Fermi energy and, 
therefore, possess the most d electrons out of the family 
of four Weyl semimetal compounds. 


IV. CHARGE CARRIERS 

In this section we investigate the properties of the 
Fermi surface. We first show the Fermi velocities of the 
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TABLE III. Mean velocities at Weyl nodes along x, y, and z 
directions in units of 10 5 m-s _1 . c and v are for the conduction 
band and the valence band, respectively. 





TaAs TaP NbAs 

NbP 

Wi 

v x 

c 

2.5 

3.1 

2.5 

3.7 


Vx 

V 

-5.2 

-5.7 

-4.8 

-5.7 


Vy 

c 

1.2 

1.5 

1.2 

1.5 


Vy 

V 

-3.2 

-3.6 

-2.0 

-3.0 


V z 

c 

0.2 

0.2 

0.1 

0.0(3) 


V z 

V 

-0.2 

-0.2 

-0.1 

-0.0(3) 

w 2 

V x 

c 

2.4 

2.3 

2.4 

2.1 


Vx 

V 

-4.3 

-4.1 

-3.3 

-3.2 


Vy 

c 

3.5 

3.0 

2.3 

2.1 


Vy 

V 

-1.7 

-2.0 

-1.2 

-1.6 


V z 

c 

4.3 

4.4 

3.7 

3.8 


V z 

V 

-1.6 

-1.5 

-1.1 

-1.0 


Weyl fermions and the band dispersion between the ener¬ 
gies of Weyl nodes and the Fermi energy for W\ and W 2 . 
Because of symmetry, we will show results for W\ and W 2 
at positive k x , listed in Table an Then the detailed stud¬ 
ies of all the electron- and hole-like Fermi surfaces found 
in the four Weyl semimetal compounds are discussed. 


A. Fermi velocity of Weyl fermion 

In Weyl semimetals, the energies of Weyl nodes are 
close to the Fermi energy. Therefore, the physical prop¬ 
erties of Weyl fermions are easy to experimentally access 
by slightly electron or hole doping the systems. The ve¬ 
locity of electron can be calculated by (1 /h)dE/dk. The 
details of how the slope of band dispersion varies be¬ 
tween the W\ or W 2 and the Fermi energy can be found 
in Fig.[8]and Fig.[9j respectively. The velocity at W\ and 
W 2 along the k x , k y , and k z directions are listed in Ta¬ 
ble HD As long as k is close enough to the Weyl nodes, 
the three-dimensional linear dispersive feature of Weyl 
cones is expected. It can be found that the velocities in 
Table HID are highly anisotropic, and so are the band dis¬ 
persions shown in Fig. [8] and Fig. [9] Consequently, the 
Fermi surface should exhibit such anisotropy. A detailed 
discussion of the Fermi surfaces will be addressed in the 
subsequent subsection. 


B. Fermi surface 

The properties of Fermi surfaces discussed below were 
calculated using discrete grids in k space. The calcula¬ 
tions were performed using a conventional unit cell. Each 
length of the grid was 0.0015, 0.0015, and 0.005 recipro¬ 
cal lattice vector along the k x , k y , and k z direction, re¬ 
spectively. Information regarding the volume and cross- 
sectional area of Fermi surfaces and the Fermi velocities 
for each compound are listed in Tabel UVl 




- 0.02 - 0.01 0 0.01 0 . 02 - 0.02 - 0.01 0 0.01 0.02 




-0.06-0.03 0 0.03 0.06-0.06-0.03 0 0.03 0.06 




FIG. 8. (Color online). Band dispersion around Wi along (a) 
k x , (b) k y , and (c) k z directions of TaAs, TaP, NbAs, and 
NbP. The k -space coordinate of each Weyl node is set to 0. 


1. TaAs 


The Fermi surfaces of TaAs are formed by both elec¬ 
tron and hole pockets, which are shown in Fig. [TOJ Ac¬ 
cording to Table HD both W\ and W 2 are below the Fermi 
energy, with the latter being closer to it. As a result, an 
electron pocket enclosing W\ or W 2 Weyl node is ex¬ 
pected. The Fermi surface of W\ is much bigger than 
that of W 2 and shows high anisotropy along k z . Some 
distance away from the end of the Fermi surface of W\, 
a hole pocket centered on the mirror plane appears and 
extends to the neighborhood of the Fermi surface of W 2 . 
Since each of the Swi and Sw2 Fermi surfaces encloses 
one Weyl node, the Fermi surface is endowed with a ±1 
Chern number^ Furthermore, the hole pocket does not 
surround any Weyl nodes, which is an outcome of the 
modulation of bands, and, consequently, carries a zero 
Chern number. The shape of the Fermi surfaces is rem¬ 
iniscent of the nodal rings in the absence of spin-orbit 
coupling. We note that different from the other com¬ 
pounds discussed below, the Fermi surfaces on two sides 
of a mirror plane are isolated instead of getting immersed, 
which is due to the distances between Weyl nodes and 
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-0.04-0.02 0 0.02 0.04-0.04-0.02 0 0.02 0.04 

k (A *) k (A *) 



FIG. 10. (Color online), (a) Top view and (b) bird’s eye view 
of TaAs Fermi surfaces in the first Brillouin zone of conven¬ 
tional unit cell with (c) the more detailed plots. Electron-like 
and hole-like Fermi surfaces are colored pink and green, re¬ 
spectively. Swi (Sw 2 ) is the Fermi surface enclosing the Wi 
(W 2 ) Weyl node(s), and S e (Sh) is an electron (hole) pocket 
without enclosing any Weyl node. 


FIG. 9. (Color online). Band dispersion around W 2 along (a) 
k x , (b) k y , and (c) k z directions of TaAs, TaP, NbAs, and 
NbP. The k- space coordinate of each Weyl node is set to 0. 


their energies relative to the Fermi energy. Furthermore, 
an electron pocket around (0.2425, 0.2425, 0) is found, 
which does not appear in other compounds. The details 
of the shapes and carrier concentrations can be found in 
Fig. Hm and Table El The maximal areas of the cross 
sections for the Fermi surfaces, which can be specifically 
measured by cyclotron experiments, are also listed in Ta¬ 
ble nvi 


2. TaP 

The Fermi surfaces of TaP are presented in Fig. [lTJ 
Contrary to TaAs, in which both energies of W\ and W 2 
are below the Fermi energy, the energy of W\ and W 2 
in TaP are below and above the Fermi energy, respec¬ 
tively shown in Table [TTJ The W 2 in TaP is enclosed by 
a hole pocket Fermi surface. Due to higher energies of 
Weyl nodes away from the Fermi energy, the volume of 
Fermi surfaces are much larger (about 10 ~ 20 times for 
Sw 1 ) than those in TaAs. Furthermore, a Fermi sur¬ 
face encloses more than one Weyl node. An electron-like 
Fermi surface encompasses a mirror pair of W\ points, 


(a) TaP (c) 




FIG. 11. (Color online), (a) Top view and (b) bird’s eye view 
of TaP Fermi surfaces in the first Brillouin zone of conven¬ 
tional unit cell with (c) the more detailed plots. Electron-like 
and hole-like Fermi surfaces are colored pink and green, re¬ 
spectively. Sw 1 (*SVv 2 ) is the Fermi surface enclosing the W\ 
(W 2 ) Weyl node(s). 
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FIG. 12. (Color online), (a) Top view and (b) bird’s eye view 
of NbAs Fermi surfaces in the first Brillouin zone of conven¬ 
tional unit cell with (c) the more detailed plots. Electron-like 
and hole-like Fermi surfaces are colored pink and green, re¬ 
spectively. Sw 1 (Sw 2 ) is the Fermi surface enclosing and 
evolved from the Wi (W 2 ) Weyl node(s), and S e (Sh) is an 
electron (hole) pocket not evolved from the Weyl node. 



FIG. 13. (Color online), (a) Top view and (b) bird’s eye view 
of NbP Fermi surfaces in the first Brillouin zone of conven¬ 
tional unit cell with (c) the more detailed plots. Electron-like 
and hole-like Fermi surfaces are colored pink and green, re¬ 
spectively. Sw 1 (Sw 2 ) is the Fermi surface enclosing and 
evolved from the Wi (W 2 ) Weyl node(s), and S e (Sh) is an 
electron (hole) pocket not evolved from the Weyl node. 


while a hole-like Fermi surface covers two mirror pairs 
of W 2 points. Compared to that in TaAs, both electron 
and hole pockets take the shape of a crescents that are 
distributed more along the k z direction. Except these 
electron and hole pockets there is no other piece of Fermi 
surfaces. The carrier concentrations and maximal areas 
of cross sections are listed in Table nvi 


3. NbAs 

The Fermi surfaces of NbAs are presented in Fig. |T21 
In NbAs, the energy of W\ and W 2 are below and above 
the Fermi energy, respectively. Three types of hole-like 
Fermi surface are found in NbAs. The first kind encloses 
four W 2 points and passes through the k z — 0 plane, 
like those found in TaP. The second kind is similar to 
those crescents elongated along the k x or k y , direction 
in TaAs. The third kind cannot be seen in Fig. [12] be¬ 
cause the Fermi energy cuts through two spin split va¬ 
lence bands at k ~(0, 0.2152, 0). Consequently, the inner 
surface is hidden inside the outer surface. In comparison 
with the W\ points in TaAs and TaP, the two neighbor¬ 
ing Wi points are closer to each other. As a result, only 
one big crescent-shaped electron pocket containing two 
W\ points is observed. Similar to the scenario of the two 
concentric-like hole pockets, there is also another elec¬ 
tron pocket hidden inside the bigger Fermi surface near 
the W\. The carrier concentrations and maximal areas 
of cross sections for each Fermi surface are listed in Ta¬ 
ble IE 


4 . NbP 

The Fermi surfaces of NbP are presented in Fig. [13] 
Similar to TaP and NbAs, the energy of W\ and W 2 in 
NbAs are below and above the Fermi energy, respectively. 
The two concentric-like hole Fermi surfaces split by spin- 
orbit coupling are also present in NbP as intersecting the 
the k z =0 plane at k «(0, 0.2042, 0), giving rise to four 
crescent Fermi surfaces in the first Brillouin zone with 
each surface containing one more surface inside. Each 
outer surface encloses one pair of W 2 points. In NbP, 
one long crescent electron-like Fermi surface resembling 
the shape of half a ring can be found. The surface be¬ 
comes thinner while deviating from the W\ and the de¬ 
tailed shape is easier to be observed in Fig. fl3lc). There 
are two additional electron-like Fermi surface inside this 
long crescent surface, S e l and S e 2. Referring to the band 
dispersion in Fig. [8] we have to carefully identify the con¬ 
centric electron pockets. Although S e 1 encloses one WI 
point, the associated band does not evolve from the band 
of Weyl nodes and, therefore, is a trivial Fermi surface. 
The size of the surface S e 2 is small and touches the limit 
for the chosen grid size. The carrier concentration and 
maximal areas of cross sections for each Fermi surface is 
listed in Table [IV] 
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TABLE IV. Electron (e) and hole (h) carrier concentrations 
(n), maximal areas (A) of cross sections on the yz, xz, and xy 
planes, and root-mean-square velocities (€>) along the k x , k y , 
and k z directions of Fermi surfaces (c/. Fig. HOD . Swi{Sw 2 ) 
is the Fermi surface enclosing and evolved from the Wi(W 2 ) 
Weyl node(s), and S e (Sh) is an electron (hole) pocket Fermi 
surface not evolved from any Weyl node. The units of carrier 
concentration, area, and velocity are 10 17 cm -3 , 10 -3 A , 
and 10 1 2 3 4 5 6 7 m • s -1 , respectively. 


TaAs 


n 

Ay Z 

A xz 

A-xy 

Vx 

Vy 

V z 


Swi(e] 

) 3.336 

6.060 

2.557 

0.639 

2.920 

1.532 

0.453 


Sw2(e] 

) 0.134 

0.244 

0.237 

0.226 

2.055 

1.579 

1.763 


Se 

0.364 

0.406 

0.406 

0.624 

0.365 

0.365 

0.412 


Sh 

2.201 

7.332 

0.310 

1.639 

4.051 

0.670 

1.372 

TaP 


n 

Ayz 

A xz 

A X y 

Vx 

Vy 

V z 


Swi(e] 

) 48.98 

23.67 

15.87 

4.687 

3.270 

1.594 

0.767 


Sw2(h) 53.17 

36.96 

13.94 

3.404 

3.736 

2.542 

1.879 

NbAs 


n 

Ayz 

A X z 

A X y 

V x 

Vy 

V z 


Swi(e] 

) 14.60 

15.97 

6.882 

1.752 

3.750 

1.370 

0.542 


Sw2(h) 29.76 

24.57 

10.29 

2.646 

4.647 

2.647 

1.765 


S e 

0.196 

1.704 

0.375 

0.060 

2.331 

1.131 

0.194 


Shi 

1.925 

3.533 

2.299 

0.470 

2.815 

1.922 

0.600 


Sh2 

0.155 

1.410 

0.037 

0.104 

3.692 

0.695 

1.887 

NbP 


n 

Ayz 

A X z 

A X y 

Vx 

Vy 

v z 


Swi(e] 

) 32.64 

31.47 

10.78 

2.805 

4.047 

1.653 

0.932 


Sw 2 (h) 82.16 

43.94 

20.15 

5.913 

4.964 

3.147 

2.551 


Sel 

4.323 

11.12 

2.928 

0.599 

4.773 

1.845 

0.501 


Se 2 

0.025 

0.219 

0.016 

0.056 

3.828 

0.466 

1.896 


Shi 

33.95 

22.59 

13.03 

3.772 

4.971 

3.294 

1.929 


V. DISCUSSION 

Several interesting properties of the charge carriers can 
be found even without the consideration of additional 
conditions, like applying magnetic or strain fields in first- 
principles calculations. We first note that the energy 
difference between W\ and W 2 is ^ 10 — 100 meV, which 
is realistic to access by doping in real crystals. The Fermi 
energy for TaAs was experimentally found^i to be located 


between W\ and W 2 (Ewi and Ew2, respectively), which 
is approximately 11.5 meV above Ew\ and, therefore, 
different from our calculations shown in Table m 

The energy difference between two Weyl nodes is di¬ 
rectly related to the minimal area of the Fermi surface 
as a function of binding energy, with a larger separation, 
\Ew 2 — E w i|, corresponding to a larger minimal Fermi 
surface. Referring to Table nn the energy differences of 
Ew 2 and Ew 1 are 0.0132, 0.0727, 0.0364, and 0.0793 
eV for TaAs, TaP, NbAs, and NbP, respectively. Recall, 
Ew 2 > Ew i- Note that while the separation in mo¬ 
mentum space of the Weyl nodes is determined by the 
strength of the spin-orbit coupling, the separation in en¬ 
ergy of the Weyl nodes can be viewed as arising from the 
dispersion of the spinless line node crossing, observed be¬ 
fore spin-orbit coupling is introduced. Amongst the four 
Weyl semimetal compounds, the relatively small value for 
TaAs explains the relatively small carrier concentration, 
e.g. 0.134xl0 17 cm -3 for Sw 2 - Similarly, the largest hole 
pocket can be found in NbP, e.g. 82.16xl0 17 cm -3 for 
Sw 2 • The rich and tunable behavior of charge carriers in 
these four semimetals strongly suggests that these com¬ 
pounds are a good playground for engineering the elec¬ 
tronic properties through doping, strain, electric fields, 
magnetic fields and/or other external perturbations. 

Furthermore, we also point out that an abrupt change 
in the Fermi surface volume may occur by tuning, for 
example, the Fermi energy. Since each pair of Weyl nodes 
are located on two sides of the mirror plane, after turning 
on the spin-orbit coupling, they are, in general, close to 
each other. Therefore, there exists a critical Fermi energy 
that can make the two Fermi surfaces, each enclosing 
the Weyl node(s) on one side of the plane, touch at one 
k point. By changing the critical Fermi energy, either 
one big surface enclosing both Weyl nodes or two small 
surfaces with each enclosing a single Weyl node can be 
formed. The Fermi surfaces of the W\ or W 2 pairs in 
TaAs is an example for the latter case while those in 
TaP, NbAs, and NbP are examples of the former case. 
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